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Abstract 

I review the potential for observing cosmic reionization using the HI 21cm line of 
neutral hydrogren. Studies include observations of the evolution of large scale struc- 
fSj ' ture of the IGM (density, excitation temperature, and neutral fraction), through 

HI 21cm emission, as well as observations of small to intermediate scale structure 
through absorption toward the first discrete radio sources. I summarize predictions 
I/"") ' for the HI signals, then consider capabilities of facilities being built, or planned, to 

jn ■ detect these signals. I also discuss the significant observational challenges. 
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1 Introduction 



The 21cm line of neutral hydrogren presents a unique probe of the evolution 
of the neutral intergalactic medium (IGM), and cosmic reionization. Furlan- 
etto & Briggs $2$) point out some of the advantages of using the HI line in 
this regard: (i) unlike Lya (ie. the Gunn- Peterson effect), the 21cm line does 
not saturate, and the IGM remains 'translucent' at large neutral fractions 
(Jl2h . And (ii) unlike CMB polarization studies, the HI line provides full three 
dimensional (3D) information on the evolution of cosmic structure, and the 
technique involves imaging the neutral IGM directly, and hence can easily dis- 
cern different reionization models (J3I: Eih. HI 21cm observations can be used 
to study evolution of cosmic structure from the linear regime at high red- 
shift (ie. density-only evolution), through the non-linear, 'messy astrophysics' 
regime associated with luminous source formation. As such, HI measurements 
are sensitive to structures ranging from very large scales down to the source 
scale set by the cosmological Jeans mass, thereby "making 21cm the richest 
of all cosmological data sets" (0). 

Early calculations of HI 21cm signal from a neutral IGM discuss a broad range 
of models for large scale HI structure, including Zeldovich 'pancakes' arising 
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in HDM models, explosive structure formation models, and CDM (|64tl34tl58l). 
Early (unsuccessful) searches for large scale structure in HI based on this 
wide range of predictions include searches for pancakes with masses > 10 14 
M at z ~ 3 (J66[ I53L Hit Hit ), and for even more massive structures (10 15 
M ) at z = 8.4 (5). Fortunately, the cosmological parameter space has been 
greatly reduced due to the advent of the concordance model of ACDM struc- 
ture formation, and this paper will assume the standard precision cosmology 
parameters (62). 

In this review I will summarize the HI 21cm probes of cosmic reionization and 
the neutral IGM. I will focus on the latest predictions for the expected signals, 
and the observational capabilities of telescopes being built, or planned, to 
detect these signals. I begin with a short basic review of the physical processes 
involved in reionization that relate directly to the HI 21cm signal. Previous 
reviews that dealt with, at least in part, the HI 21cm signal from cosmic 



reionization, include (J4| |4J| 1441 1 . and in the context of the Square Kilometer 
Array (SKA) 

See also articles by Furlanetto, Morales, and Peterson 

in this volume. 



2 The physics of the neutral IGM 



The physics and equations of radiative transfer of the HI 21cm line through 
the neutral IGM have been considered in detail by many authors fl58t |g; l22t 
SOU lei 00, and I review only the basic results here. 



In analogy to the Gunn-Peterson effect for Lya absorption by the neutral 
IGM, the optical depth, r, of the neutral hydrogen to 21cm absorption for our 
adopted values of the cosmological parameters is: 

r = J'T'trn, \ ~ 0.0074^(1 + S )(1 + (1) 



This equation shows immediately the rich physics involved in studying the HI 
21cm line during reionization, with r depending on the evolution of cosmic 
over- densities, 5 (predominantly in the linear regime), the neutral fraction, 
xhi (ie. reionization), and the HI excitation, or spin, temperature, T s . 

In the Raleigh- Jeans limit, the observed brightness temperature (relative to 
the CMB) due to the HI 21cm line at a frequency v = vq/(1 + z), where 
uq = 1420.40575 MHz, is given by: 

T B w — — rw 7(l+5)x HI (l j—^l + z) 1 mK, (2) 

L + Z Is 
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The conversion factor from brightness temperature to specific intensity, L,, is 
given by: I v = ix^xii+z))^ 8 ' or * n more common units: I v — 7.1 x 10 4 (1+2:)~ 2 Tb 
Jy sr" 1 . Equation 2 shows that for T s ~ T C mb one expects no 21cm signal. 
When T s » T C mb, the brightness temperature becomes independent of spin 
temperature. When Ts << Tcmb, we expect a strong negative (ie. absorption) 
signal against the CMB. 



An important point to keep in mind is that, for the 21cm experiments being 
considered, the signal being observed corresponds to large scale structure, not 
individual galaxies. This point is demonstrated in Fig 1, which shows the 
expected signal in Jy versus redshift for different HI masses, along with the 
expected sensitivity of current and future radio telescopes. For simplicity, I 
have assumed a single line width of 300 km s -1 , which would correspond to 
individual virialized galaxies, or to large scale structures just separating from 
the Hubble flow. Fig 1 shows that current instruments, such as the GMRT, 
can detect large galaxies (~ 10 10 M in gas) out to only modest redshifts (z < 
0.3), even in long integration times. At redshifts corresponding to reionization 
(z > 6), future large area low frequency radio telescopes, such as LOFAR, and 
eventually the SKA, will still be limited to studying large scale structure (HI 
masses > 10 11 to 10 12 M ). Fortunately, the entire IGM is made up of neutral 
hydrogen prior to reionization, and the large scale structure detected is not 
just density enhancements (ie. protoclusters), but also structure induced by 
reionization itself (the 'bubble machine' of HII regions), and possibly spatial 
variations in the spin temperature. 
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65f l show that the HI excitation temperature will equilibrate with the CMB 
on a timescale ~ year, in absence of other effects. However, collisions and 
resonant scattering of Lya photons can drive T$ to the gas kinetic temperature, 



1 + y c + y Lya 



In this equation, y c represents collisional excitation of the hyperfine transition, 
which couples T? to the gas kinetic temperature Tk ■ The coupling coefficient 
y c oc n,ff, and (|79[) shows that, for the mean IGM density, collisional coupling 
between Ts and Tk becomes significant for z > 30. The third term in equ 3 
corresponds to the Wouthuysen-Field effect, in which resonant scattering of 
Lya photons couples the spin temperature to T K ffTol : EH : H; |33| 4A; loa ). 



f)43l ) show that this latter mechanism will be important when the radiation 
background at the Lya frequency satisfies J a > 9 x 10~ 23 (1 + z) erg cm -2 s~ x 
Hz -1 sr" 1 , or about one Lya photon per every two baryons at z = 8 (@). 

The interplay between the CMB temperature, the kinetic temperature, and 
the spin temperature, coupled with radiative transfer, lead to a number of 
interesting physical regimes for the HI 2fcm signal. (Q;[j|) suggest the following 
plausible regimes: 

At z > 200 residual free electrons couple Tcmb and T K through Thompson 
scattering and subsequent gas collisions, while the density is high enough to 
equilibrate Tk and T$. In this case T$ = Tqmb and there is no 21cm signal. 

At z ~ 30 to 200, the residual ionization fraction and density is too low to 
couple T K to T(jm b 5 so the gas cools adiabatically, with temperature falling 
as (f + z) 2 , ie. faster than the (1+z) for the CMB. Hence the gas becomes 
colder than the CMB. However, the mean density in this redshift range is still 
high enough to provide some coupling between T s and T K through collisions, 
and the HI 21cm signal might be seen in absorption against the CMB. In this 
regime the HI fluctuations are still evolving linearly, essentially following the 
dark matter (3 BUI). 



At z ~ 20 to 30, the situation starts to become complex. Collisions can no 
longer couple Tk to Ts, and Ts again approaches Tcmb- However, we might 
also expect the first luminous structures (Pop III stars or mini-quasars), at 
least near the end of this redshift range. The Lya photons from these objects 
would induce local coupling of T K and To, thereby leading to some 21cm 
absorption regions. On the other hand, (yl HI) point out that these same 
pho tons, and more importantly, any Xrays from the first luminous sources 



( 141). could lead to local IGM warming above Tcmb well before reionization. 

±£ 
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Energetically it takes 13.6 eV per baryon to ionize the IGM, but only 0.005 ^ 1+2 ^ 
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eV to warm the IGM about Tcmb- Hence one might expect a 'patch work' of 
regions with no signal, absorption, and perhaps emission, in the 21cm line. 




as well as 
such that 



At z ~ 6 to 20 all the physical processes come to play. The IGM is being 
warmed by the resonant scattering of Lya photons (|43f). and penetrating, ion- 
izing hard Xrays, from the first galaxies and black holes 
by weak shocks associated with structure formation (25| 2 
T K is likely larger than T C mb globally (25). Likewise, these objects are reion- 
izing the universe, leading to a fundamental topological change in the IGM, 
from the linear evolution of large scale structure, to a bubble dominated era 
of HII regions (J26h . It is this regime that most (although not all), theoretical 
work has been exploring, and the expected HI 21cm signal is rich. Finally, 
after reionization, z < 6 to 10, the IGM is fully ionized (xhi < 10 -4 ), and the 
21cm signal is gone. 



3 HI 21cm probes of the evolution of the IGM 



3.1 Global HI signature 



In an idealized sense, if the universe reionized very rapidly everywhere, at the 
same time, one would expect a global (ie. full sky) step in the background 
temperature at the frequency corresponding to the redshifted 21cm line (60). 
However, if nothing else, cosmic variance will lead to different ionization red- 
shifts in different regions, and we now know that reionization is likely an 
extended process in time, thereby smoothing out the expected signal in space 
and time. On the other hand, since this is an all sky signal, the sensitivity of 
the experiment is independent of telescope collecting area, and the experiment 
can be done using small area telescopes at low frequency, with well controlled 
frequency response. 



The most recent modeling of the expected global HI signal is presented in ()27l ) 
(see also (J25f)). who employ the numerical simulations of (j29h. They include 
both a fast reionization model, plus more complicated, ad hoc reionization 
histories. 

Results from these calculations are reproduced in Fig 2a. The expected signal 
peaks at roughly 20 mK above the foreground. In their fiducial (fast reioniza- 
tion) model, the peak occurs at z ~ 10. They also predict a negative signal, 
meaning absorption against the CMB, at higher redshift, prior to IGM warm- 
ing, but allowing for Lya resonant scattering (ie. the era of 'Lya coupling 
(|3|)). The shaded region shows the potential system thermal noise of a well 
calibrated low frequency experiment. 
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Fig. 2. Left: Global (all sky) HI signal from reionization |27h . The shaded region 
shows the expected thermal noise in a carefully controlled experiment. Right: The 
simulated spectrum of a source with S120 = 20 mJy at z = 10 using a spectral model 
based on the powerful radio galaxy Cygnus A and assuming HI 21cm absorption by 
the IGM (^3). Thermal noise has been added using the specifications of the SKA 
and assuming 10 days integration with 1 kHz wide spectral channels. The onset of 
absorption by the neutral IGM is seen at 129 MHz, corresponding to the HI 21cm 
line at z = 10. 



(|27l ) point out that detecting this signal against the mean non-thermal fore- 
ground radiation may be difficult. The foreground is the sum of relatively 
smooth Galactic synchrotron emission and discrete distant radio galaxies, with 
a typical spectral index, a ~ —0.8, where a is defined as: S u oc v a . The fore- 
ground temperature behaves roughly as Tfg ~ 100 ( 200 ^ 1Hz )" 2,8 , in the coldest 
regions of the sky, and can be an order of magnitude higher in the Galactic 
plane. Hence, the expect change in Tb is at most ~ 10 -4 that of the mean 
foreground signal over the frequency range ~ 100 to 200 MHz. 



3.2 Large scale structure 



3.2.1 Power spectra and tomography 

The majority of theoretical studies have focused on predicting the HI fluctua- 
tions in the IGM on scales of arcmins to a few degrees both during ffFsL li^i Irifit 
3; H 3 El IS HQ , and prior to (§ 3 E3) , reionization. Most studies 



have focused on the diffuse IGM, although a few studies have also considered 
the signal expected due to the clustering of collapsed structures (|36f). 

The HI 21cm signal from the IGM during reionization can be predicted both 
analytically, using a standard Press-Schechter type of analysis of linear struc- 



6 



d 

-—I 

s 




-2 -1 1 2 
x (arcmin) 



-2-10 1 2 
8 (arcmin) 



-2-10 1 2 
s (arcmin) 



-2-10 1 
log <5T» (mK) 



Fig. 3. The simulated HI 21cm brightness distribution during reionization at z = 12, 

9, 7 Hii). 



ture formation, plus some recipes to approximate non- linear evolution ([571 : l7£ 
El IfJ H3), or through the use of numerical simulations (@; HF 



Fig 3 shows the expected evolution of the HI 21cm signal during reionization 
based on the numerical simulations of (j^). They find that the mean HI signal 
is about Tb ~ 25 mK prior to reionization, with fluctuations of only a few 
mK on arcmin scales due to linear density evolution. In this simulation, the 
HII regions caused by galaxy formation during reionization are seen in the 
redshift range z ~ 8 to 10, reaching scales up to 2' (frequency widths ~ 0.3 
MHz => 0.5 Mpc physical size). These regions have (negative) brightness 
temperatures up to 20 mK relative to the mean HI signal. This corresponds 
to 5/iJy beam -1 in a 2' beam at 140 MHz. 

We compare this signal to the sensitivity of a radio array. The point source 
rms sensitivity (dual polarization) in an image from a synthesis radio telescope 
is given by: 
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where Au is the channel width in kHz, t is the integration time in hours, T sys 
is the total system temperature, A is the collecting area of each element in 
the array, e is the aperture efficiency, and iV is the number of elements (eAN 
= total effective collecting area of the array). At low frequency the value 
of T sys is dominated by the non-thermal foreground, and again behaves as 
T sys ~ 100 ( 200MHz ) ~ 2 8 The beam FWHM is given by the Fourier transform 
of the uv coverage, assuming equal weight for each visibility. 

Consider an array with an effective collecting area of 1 square kilometer at 140 
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MHz, distributed over 4 km, a system temperature of 300 K, a channel width 
of 0.3 MHz, and integrating for one month. The rms sensitivity is then 1.3/iJy 
beam -1 , with a beam FWHM ~ 2'. This square kilometer array will be just 
adequate to perform true three dimensional imaging of the average structure 
of the IGM during reionization. 

Unfortunately, the nearer term low frequency 'path finder' arrays will have 
< 10% the collecting of the SKA, and will likely not be able to perform such 
direct 3D imaging (Table 1). However, these near term experiments should 
have enough sensitivity to perform power spectral analyses of the HI 21cm 
fluctuations. For power spectral analyses the sensitivity is greatly enhanced 
relative to direct imaging due to the fact that the universe is isotropic, and 
hence one can average the measurements in annuli in the uv-plane, ie. the 
statistics of fluctuations along a uv annulus are equivalent. The width of in- 
dependent uv annuli is set by the primary beam diameter, D, or Al = 27i-^-. 
The situation is analogous to COBE and WMAP. COBE lacked the sensitiv- 
ity to make a true image of the CMB fluctuations, but was able to make a 
robust determination of the power spectrum of fluctuations. WMAP was able 
to produce a proper image of the CMB. 

Many authors have considered the power spectrum of brightness temperature 
fluctuations in the HI 21cm line (0; ES B B H H 0) • It can be shown 
that the noise power spectrum, in mK 2 , in standard spherical harmonic units, 
;a C/^) for an array with uniform coverage of the uv plane, is given by: 



AutfX 



■max 



where f c = areal covering factor of the array = N ant ^^-, and A tot is the 
total area of the area defined by the longest baseline. For example, for the 
logarithmic antenna spacing of the VLA in the smallest configuration, the 
VLA-VHF system has f c = 0.017 on baselines out to 1 km, and f c = 0.17 
out to 0.25 km, for a full synthesis observation. For comparison, the SKA will 
have 30% of it's collecting area inside 1 km, or f c = 0.3 out to 1 km. 

Most predictions of the HI 21cm power spectrum during reionization published 
thus far parallel the 2D calculations applied to the CMB. Ultimately, higher 
sensitivity, and physical insight, will come through full 3D analysis of the 



HI 21cm fluctuations during reionization (48). Of course, at the smallest I, 
approaching the primary beam size Al, the signal-to- noise is limited by cosmic 
variance to essentially the number of synthesized beams per primary beam. 

As an example, Fig 4 shows the predicted evolution of the HI 21cm power 
spectrum using an analytic calculation of structure formation, including a 
model of (uncorrelated) HII region formation during reionization (25[ 7^) (see 
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Fig. 4. Predicted HI 21cm brightness temperature power spectrum during reioniza- 
tion at z = 10 ({22I : F^ ). The green dotted line is the model prediction including 
only linear density fluctuations. The black line solid includes the effect of HII re- 
gions during reionization. The red long dash line shows the sensitivity of pathfinder 
experiments like LOFAR, while the blue short dash curve shows the sensitivity for 
the SKA, assuming log / bins. The cutoff at low / is set by the primary beam. 



also (p; |^2t IfiSD f° r analytic solutions, and (j25t for predictions based on 
numerical simulations). Fig 4 includes the power spectrum arising from linear 
density fluctuations following the dark matter, plus the effect of reionization. 
The signature of reionization can be seen as a bump in the power spectrum 
above the density-only curve due to the formation of HII regions. In this cal- 
culation, the rms fluctuations at z — 10 peak at about 10 111K rms (observing 
frequency = 130 MHz) on scales / ~ 3000 1 . Note that in this calculation the 
neutral fraction as a function of redshift is a free parameter, and is assumed 
to be 0.5 at z = 10. (0) point out that there is no very large scale information 
contained in the HI 21cm fluctuations, and that power spectral analyses are 
best done on scales ranging from / = 100 to a few thousand, ie. on angular 
scales < few degrees. 



Also included are the noise power spectra for near and longer term radio 
arrays. Typical noise values for near-term arrays are ~ 1 to 10 mK rms in 
the range I = 10 3 to 10 4 , for long integrations. This should be adequate to 
determine the HI 21cm power spectrum during reionization. The SKA will go 
an order of magnitude deeper, providing an accurate measure of the HI power 
spectrum, and its evolution with redshift. 



1 or u —— 7t — —— 500, or 1 

Z7T 



~ ^r- = 3.6', or baselines, b 



\u = 1.2 km, or wavenumber 



/ 



Mpc 1 = 0.3Mpc 1 (comoving). 
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(J2 E|) have shown that gas peculiar velocities, ie. infall into superclusters 
and along filaments (analogous to the Kaiser effect for galaxy clustering), 
increases the expected rms fluctuations by about a factor two. Also, clustering 
of luminous sources (ie. biased galaxy formation) implies a characteristic scale 
for the HII regions of 1 to 10 Mp c (comoving), and leads to another factor two 



in fluctuation strength (26; l57f ) over the random distribution of HII regions 
assumed by (|77l). Specifically, this later effect moves the fluctuation power to 
larger scales. 

(0; IE 0) consider the effect of spectral channel width on the power spectral 
analysis. Bandwidth enters via the finite thickness of the HII regions and the 
density fluctuations. They conclude that a bandwidth of about 0.2 MHz is 
optimal for power spectral studies, although the rms of the signal changes by 
only a factor of about two going from 0.1 to 0.5 MHz channel width. 



f)25l ) show how the power spectra change with different reionization models, 
including outside-in versus inside-out models, uniform versus HII region dom- 
inated reionization, and double reionization. They point out that HI 21cm 
power spectral studies should provide powerful constraints on the process of 
reionization, with different reionization scenarios clearly discernable in their 
power spectral signature. Note that this is unlike CMB large scale polarization 
measurements, which provide reasonable constraints on the total Thompson 
scattering optical depth through the universe back to recombination, but only 
marginal constraints on the reionization history (J35l). 



Lastly, (|6j 1491 ) point out that if full 3D tomography of the evolution of IGM 
structures could be made, then one might also constrain the geometry of uni- 
verse through the standard Alcock-Pacinsky effect, ie. the difference in struc- 
ture evolution observed in angle versus along the line-of-sight, due to the 
non-Euclidean geometry of the univese. 



3. 2. 2 Clustering of minihalos 

A point of debate has been the fraction of HI in collapsed objects, as opposed 
to the diffuse IGM (|36h . This fraction has a complex depence on structure 
formation history, and (52; 28) conclude that the majority of the HI during 
reionization will remain in the diffuse phase. However, at times when T$ ap- 
proaches T C mb in the diffuse IGM, one might expect HI 21cm emission due to 
clustering of minihalos 2 , in which T$ > Tqmb simply due to the virialization 

2 Minihalos have masses < 10 7 M©, and virial temperatures < 10 4 K. These halos 
cannot cool via atomic hydrogen lines, and therefore 'saturate' at over-densities of 
5 ~ 100. They will not form stars unless an alternate cooling mechanism can be 
found. One candidate is molecular hydrogen, which is an effective coolant down to 
virial temperatures of 100 K (0). 
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process (25). 



(|36l ) have considered the HI 21cm power spectrum due to clustering of miniha- 
los at high redshift. For a beam size of 1', and channel width of 0.2 MHz, they 
predict 3a brightness temperature fluctuations due to clustering of minihalos 
~ 7mK at z = 8.5, decreasing to 2mK at z = 20. 



3.2.3 Other tests 

(jl6h discuss circular polarization of the 21cm line due to zeeman splitting 
during reionization. Such polarization would be detectable with the SKA if the 
IGM field strength were 100 fiG. ( 
HI 21cm fluctuations and sub-degree scale CMB secondary anisotropies due 
to Thompson scattering in ionized regions during reionization at / > 1000 
(kinetic SZ effect). A natural anticorrelation would be expected between the 
neutral and ionized regions. And lastly, (J6l|) discuss delensing of the CMB 
using the HI 21cm fluctuations to pin-point density inhomogenaeities along 
the line of sight (z ~ 10 to 200). This could prove important for decoupling 
the E and B mode mixing that occurs due to lensing, possibly recovering the 
intrinsic (inflationary) B polarization. 



mi 



consider the anti-correlation between 



3.2.4 Cosmic Stromgren spheres 

While direct detection of the typical structure of HI and HII regions may be 
out of reach of the near-term 21cm telescopes, there is a chance that even this 
first generation of telescopes will be able to detect the rare, very large scale 
HII regions associated with luminous quasars near the end of reionization. 

The presence of cosmic Stromgren spheres around the highest redshift SDSS 
QSOs has been deduced from the observed difference between the redshift 
of the onset of the GP effect and the systemic redshift of the host galaxy 



f)69l ; l72l; l7ll; 1671) (although cf. (|50j)). The physical size of these spheres is ~ 5 
Mpc, or an order of magnitude larger than the typical spheres expected from 
clustered galaxy formation, due to the extreme luminosity of the sources (~ 
10 14 L ). If the neutral fraction remains substantial to relatively low redshift, 



f(HI) > 0.1 (|46t l72t 1711 ). then it is plausible to search for these regions as 
'holes in the sky', or regions of negative brightness in the redshift 21cm line 
frequency. The expected signal is ~ 20mK x f(HI) on a scale ~ 10' to 15', 
with a line width of ~ 1 to 2 MHz |73h. This corresponds to 0.5 x f(HI) mJy 
beam -1 , for a 15' beam. 

Searching for HI signals around known high redshift QSOs has a number of 
major observational advantages: (i) the exact location in frequecy, and RA and 
DEC, are well determined, thereby limiting the search space dramatically, and 
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RIGHT ASCENSION [J2000] 




Fig. 5. Simulated image and spectrum of a Stromgren sphere surrounding a lu- 
minous SDSS QSO at the end of reionization for a 250 hour integration with the 
VLA-VHF system (0). The contour levels are: -0.45, -0.3, -0.15, 0.15, 0.3, 0.45 mJy 
beam -1 , and the negative signal has been inverted to appear positive. 

(ii) there is already evidence for the existence of the features based on the Lya 
spectra. The VLA-VHF system has been specifically design to search for these 
spheres around known SDSS quasars at z ~ 6 to 6.4 (j86l ). Fig 5 shows the 
simulated image and spectrum of such a Stromgren sphere around an SDSS 
QSO assuming f(HI) = 1, as observed with the VLA-VHF system. While it is 
unlikely the mean neutral fraction is this high at such a low redshift, in 250 
hours the VLA-VHF system should set the first direct constraints on neutral 
fraction of the IGM, at the level f(HI) > 0.2, as well as easily rule out more 
extreme models, such as T s << T CMB . (J74|) present similar simulations for 
the MWA and the SKA. 



f)39f l also discuss spectral dips due to large HII regions around luminous 
quasars (> 2 x 10 10 L Q ) during the EoR. Their Fig 2 shows a spectrum in 
a typical 10' beam from 100 to 180 MHz. They predict on average one rel- 
atively deep (-2 to -4 mK) dip per LoS on this scale. They emphasize that 
these spectral features may be easier to detect than spatial fluctuations of a 
similar magnitude, due to the fact that spatial confusion is highly structured 
on arcsec to arcmin scales, while the spectral confusion should be smooth over 
10's of MHz. 



(J74]) perform a similar calculation using the evolution of the bright QSO lumi- 
nosity function to predict the number of HII regions around active QSOs at 
z > 6. They conclude that there should be roughly one SDSS-type HII region 
around an active QSO (physical radius > 4Mpc) per 400 deg 2 field per 16 Mhz 
bandwidth at z ~ 6, and one R > 2 Mpc region at z ~ 8. They also point 
out that the recombination time may be longer than the Hubble time, so that 
fossil HII regions may be observed around non-active AGN. Assuming a duty 
cycle ~ ^^TX ~ ?S ~ °- 01 > where < > is the fiducial lifetime of QSO 
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activity (75), and tn(z) is the age of the universe, leads to a factor 100 more 
fossil HII regions per FoV. 

3.2.5 Beyond reionization 

A number of authors have considered the brightness temperature fluctuations 
due to the 21cm line from the neutral IGM at redshifts prior to reionization, 
z > 20. (0; 0) predict the power spectrum of fluctuations during the era when 
the Lya photons from the first luminous objects couple T s to T K locally, but 
the universe remains substantially neutral [z ~ 20 to 30), and Tk < Tomb in 
the diffuse IGM. Brightness temperature fluctuations can be due to emission 
from clustered minihalos, plus enhanced by absorption against the CMB by 
the diffuse IGM, with an rms ~ lOmK for I ~ 10 5 , due to a combination of 
linear density fluctuations, plus poisson ('shot') noise and biasing in the Lya 
source (ie. galaxy) distribution. 

(Jiil ) go even further in redshift, to z > 50 to 200. In this regime the HI 
generally follows linear density fluctuations, and hence the experiments are 
as 'clean' as CMB studies, and T K < T C mb, so a relatively strong absorption 
signal might be expected. (42) also point out that Silk damping, or photon 
diffusion, erases structures on scales I > 2000 in the CMB at recombination, 
corresponding to comoving scales = 22 Mpc. The HI 21cm measurements can 
explore this physical regime at z ~ 50 to 300. The predicted rms fluctuations 
are 1 to 10 mK on scales I = 10 3 to 10 6 (0.2° to 1"). These authors point out 
that, due to sensitivity to large I, and the 3D nature of the information, HI 
21cm power spectral studies in this epoch "contain an amount of information 
orders of magnitude larger than any other cosmological probe". These data 
could provide the best tests of non-Gaussianity of density fluctuations, and 
for constrainting the tilt, or running power law index, of mass fluctutions to 
large I, providing important tests of inflationary structure formation. (jHif l also 
suggests that a large global signal, up to -0.1 K, might be expected for this 
redshift range. 

Two experiments are being planned to cover this low frequency range, the long 
wavelength array (LWA) in New Mexico, and in the longer term, LUDAR is 
being considered for the back side of the moon. 

3.3 Absorption toward discrete radio sources 

Observing HI 21cm emission from the EoR implies studying large scale struc- 
ture (cluster scales and larger). A number of groups have recently considered 
the possibility of studying smaller scale structure in the neutral IGM by look- 
ing for HI 21cm absorption toward the first radio-loud objects (AGN, star 
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forming galaxies, GRBs) (fib!). 



(0) use numerical simulations to predict the HI 21cm absorption profile of 
the 'cosmic web' prior to reionization. For example, for a source at z — 10, 
they predict an average optical depth due to 21cm absorption of about 1%, 
corresponding to the 'radio Gunn-Peterson effect'. They also find about five 
narrow (few km/s) absorption lines per MHz with optical depths of a few to 
10% (Fig 2b). These latter lines are equivalent to the Ly a forest seen after 
reionization, and correspond to over-densities evolving in the linear regime 
(5 < 10). While significant questions remain about simulating the thermal 
state of the IGM to such detail, the simple point remains that, while the 
Lya lines from such structures are highly saturated, the (much) lower 21cm 
oscillator strength makes the IGM translucent prior to reionization. 



(|23l ; |5l|) predict a similar HI 21cm absorption line density due to gas in mini- 
halos as that expected for the 21cm forest. (0) shows that the presence or 
absence of HI absorption by minihalos, or the cosmic web, would be a telling 
diagnostic of early IGM heating mechanisms. (23) also consider the expected 
21cm absorption profiles for proto-disk galaxies. While such absorption lines 
will be rare (10 4 times less frequent than the 21cm forest lines), the optical 
depths may be large enough that the lines could be observed toward faint ra- 
dio sources, in particular, gamma ray burst radio after-glows within the host 
galaxy (|37fl. 

An important caveat in these calculations is the assumption of radio loud 
sources during the EoR. This question has been considered in detail by |l2l 
libl liih. They show that current models of radio-loud AGN evolution predict 
between 0.05 and 1 radio sources per square degree at z > 6 with Sisomhz > 6 
mJy, adequate for EoR HI 21cm absorption studies with the SKA. 



4 Observational challenges 



4-1 Foregrounds 



It has long been recognized that the HI 21cm signal from reionization must be 
detected on top of a much larger non-thermal (synchrotron) signal from fore- 
ground emission. This foreground includes discrete radio galaxies, and large 
scale emission from our own Galaxy. The expected HI signal is about 10~ 4 of 
the foreground emission. 

Many groups have considered the effects of the foregrounds on HI 21cm EoR 
imaging and power spectral studies. (20) show that, even if point sources 
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can be removed to the level of 1 /iJy, the rms fTuctutions on spatial scales 
< 10' (I > 1000) due to residual radio point sources will be > 10 mK just 
do to Poisson noise, increasing by a factor 100 if the sources are strongly 
clustered. Clearly this conclusion depends on the extrapolation of radio source 
populations to < 1/xJy, but under reasonable assumptions, even in the Poisson 
case the residual source fluctuations will be comparable to the HI signal. 

This calculation has led many groups to consider removal of foregrounds in 
the spectral domain. The important point is that the foregrounds should be 
relatively smooth in frequency, predominantly the sum of power-law, low fre- 
quency non-thermal spectra, or perhaps gently curving on spectral scales 
> 10's MHz. The HI signal should show significant structure on sub-MHz 
scales, corresponding to the typical size scale for features during reionization. 



A number of complimentary approaches have been presented for foreground 
removal. (27; 6^) consider removal of the foregrounds through fitting of smooth 
spectral models (power laws or low order polynomials in log space) to the 
observed visibilities or images. (|46[ liih present a 3D Fourier analysis of the 
measured visibilities, where the third dimension is frequency. They show that 
the different symmetries in this 3D space for the signal arising from the noise- 
like HI emission, versus the smooth (in frequency) foreground emission, can be 
a powerful means of differentiating between foreground emission and the EoR 
line signal. (|57l |(J l76h perform a similar analysis, only in the complementary 
Fourier space, meaning cross correlation of spectral channels. They show that 
the 21cm signal will effectively decorrelate for channel separations > 1 MHz, 
while the foregrounds do not, and this too can be a means of separating the 
foregrounds from the desired 21cm fluctuations. The overall conclusion of these 
methods is that spectral decomposition should be adequate to separate non- 
thermal foregrounds from the HI 21cm signal from reionization at the mK 
level. 



However, there are two potentially more insideous observational challenges 
relating to non-thermal foregrounds beyond simple in-beam confusion. First, 
there can be frequency dependent sidelobes from confusing sources in the very 
wide fields being considered. A particularly problematic effect has been the 
recent observation that telescope spectral response (ie. bandpass) may depend 
on position of a given source in the primary beam (j54l ). The origin of this effect 
remains uncertain (perhaps relating to scattering and blocking structures), but 
it will lead to frequency dependent sidelobes that potentially affect widefields 
and are not removed through normal (on-axis) bandpass calibration, and may 
require a spatially dependent bandpass calibration. 

And second is polarized structure in the diffuse emission in the field. Wide field 
polarization has been seen with low frequency observations with the WSRT at 
330 MHz and below, and is thought to arise due to differential Faraday depth 
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through the ISM (|31l). Hence, great care is required to control, and calibrate, 
the wide field polarization response of the primary elements. 



4-2 Ionosphere 



A second potential challenge to low frequency imaging over wide fields is phase 
fluctuations caused by the ionosphere. These fluctuations are due to index of 
refraction fluctuations in the ionized plasma, and behave as A0 oc v~ 2 . Mor- 
ever, the typical 'isoplanatic patch', or angle over which a single phase error 
applies, is a few to 10 degre es ( physical scales of 10 's km in the ionosphere), 
depending on frequency (jl7t 1400 . Fields larger than the isoplanatic patch will 
have multiple phase errors across the field, and hence cannot be corrected 
through standard (ie. single solution) phase self-calibration techniques. 



New wide field self-calibration techniques, involvin g multip le phase solutions 
over the field, or a 'rubber screen' phase model (|l7t I32T ) . are being devel- 
oped that should allow for self-calibration over wide fields. However, even 
self-calibration techniques will be insufficient to overcome the very rapid vari- 
ations caused by the occasional ionospheric storm, or traveling ionospheric 
disturbances. 

At very low frequency (< a few MHz), one approaches the plasma frequency 
of the ionosphere, and the optical depth increases dramatically, precluding 
observations from the ground. 



4.3 Interference 



Perhaps the most difficult problem facing low frequency radio astronomy is 
terrestrial (man-made) interference. The relevant frequency range corresponds 
to 7 to 200 MHz [z = 200 to 6). These are not protected frequency bands, and 
commercial allocations include everything from broadcast radio and television, 
to fixed and mobile communications. 

Many groups are pursuing methods for RFI mitigation and excision (see ex- 
tensive references at fjiol )). These include: (i) using a reference horn, or one 
beam of a phased array, for constant monitoring of known, strong, RFI sig- 
nals, (ii) conversely, arranging interferometric phases to produce a null at the 
position of the RFI source, and (iii) real-time RFI excision using advanced 
filtering techniques in time and frequency, of digitized signals both pre- and 
post-correlation. This requires very high dynamic range (many bit sampling), 
and very high frequency resolution. 
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One obvious advantage of HI absorption experiments (section 3.3) is that they 
can be done using long baselines (10's to 100's of km). Long baselines lead 
to decorrelation of the terrestrial signal due to fringe tracking of the celestial 
source. However, the fringe rates on the short baselines required for the HI 
emission experiments (< few km) are low, such that decorrelation due to fringe 
winding will not be a very effective RFI filter on the shorter baseline. 

In the end, the most effective means of reducing interference is to go to the 



remotest sites. (|85t 1811 ) have selected sites in remote regions of Western Aus- 
tralia, and China, respectively, because of known low RFI environments. Of 
course, the ultimate location would be the back side of the moon. 



5 Telescopes 



Many programs have been initiated to study the HI 21cm signal from cosmic 
reionization, and beyond. These are summarized in Table 1. A number of 
different approaches are being taken, both techinically and in order to address 
different aspects of the problem (jiih. 

The largest near-term efforts are the Mileura Wide Field Array (MWA), the 
Primeval Structure Telescope (PAST), and the Low Frequency Array (LO- 
FAR). These are being optimized to study the power spectrum of the HI 
21cm fluctuations, although in principle they will be able to image the larger 
HII regions during reionization. The VLA-VHF system is designed specifically 
to set limits on the HII regions around z ~ 6 to 6.4 SDSS QSOs, although 
it should also constrain the late-time power spectrum. In the long term the 
Square Kilometer Array should have the sensitivity to perform true three di- 
mensional imaging of the neutral IGM in the 21cm line during reionization. 
And at the lowest frequencies (< 50 MHz), the Long Wavelength Array (LWA), 
and eventually the lunar array (LUDAR (jH; are being designed for the 
higher z signal, prior to reionization. 

As discussed above, the technical challenges are many. Use of spectral decom- 
position to remove the foregrounds requires careful control of the synthesize 
beam as a function of frequency, with the optimal (although difficult) solution 
being a telescope design where the synthesized beam is invariant as a function 
of frequency ([901) . High dynamic range front ends are required to avoid satu- 
ration in cases of strong interference, while fine spectral sampling is required 
to avoid Gibbs ringing in the spectral response. The polarization response 
must be stable and well calibrated. Calibation in the presence of a significant 
ionospheric phase screen requires new wide field calibration techniques. The 
very high data rate expected for many element (> 10 3 ) arrays requires new 
methods for data transmission, cross correlation, and storage (85). At the low- 
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Table 1 

HI 21cm Experiments 



Preq Area B maa; Site Type FoV Date Goal Ref 

MHz 10 4 m 2 km deg 



GMRT 


150 


- 165 


3.7 


10 


India 


Parabola 


4 


2000 


StromSph 


(89) 


PAST 


50 


- 200 


7 


2 


China 


Dipolc 


10 


2006 


PowSpec 


w 


VLA-VHF 


180 


- 200 


1.3 


1 


USA 


Parabola 


1 


2006 


StromSph 


(86J 


Mark I 


100 


- 200 







Aus 


Spiral 


180 


2006 


Global 


(W) 


MWA-LFD 1 


80 


- 300 


3 


1.5 


Aus 


Dipole 


20 


2007 


PowSpec 


(85J 


LOFAR 2 


115 


- 240 


10 


2km=40% 


NL 


Tiles 


20 


2007 


PowSpec 


(82) 


LWA 


10 


- 88 


10 


5km=30% 


USA 


Dipole 


20 


2008 


Pre-reion 


m 


SKA 


100 


- 200 


100 


5km=50% 


? 


Dipole 


20 


2015 


Imaging 




LUDAR 


5 


- 50 


1000 


100 


Moon 


Dipole 


60 


2020 


pre-reion 





An order of magnitude increase in area is planned after testing the Low Frequency Demonstrator. 
2 LOFAR will also have a low frequency component between 30 and 80 MHz. 



est frequencies, < 20MHz or so, where we hope to study the pre-reionization 
IGM, phase fluctuations and the opacity of the ionosphere becomes problem- 
atic, leading to the proposed LUDAR project on the far side of the moon. 
The far side of the moon is also the best location in order to completely avoid 
terrestrial interference. 
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